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Abstract 

The completely general radiative corrections to lowest order, including the final and initial state 
radiations, are studied in proton- antiproton annihilation to an electron-positron pair. Numerical 
estimates have been made in a realistic configuraton of the PANDA (FAIR) detector for the proton 
time-like form factors measurements. 



* ong@ipno.in2p3.£r 



1 



I. INTRODUCTION 



Precise polarization measurements of the proton electromagnetic form factor [1] confirms 
the dependence up to Q"^ = 8.5 GcV^, of the ratio R = ^Ge/Gm (where /x is the proton's 
magnetic moment and Ge and Gm arc the electric and magnetic proton form factors) show- 
ing an approximately linear decrease of R with Q"^. This fact is in disagreement with results 
obtained from a new Rosenbluth cross section measurement [2] and suggests that the source 
of the discrepancy is not simply experimental. Recently, there has been a revival of interest 
in this subject [3]. Some theoretical works [4-6] have investigated the two- photon exchange 
corrections to the lowest order QED. This effect has been shown to resolve partially the 
discrepancy [3,7]. It is well established that the Rosenbluth method is much more sensitive 
to the radiative corrections than the polarization method. Until now, intense theoretical 
activities to evaluate the radiative corrections to elastic electron-proton scattering, which 
include higher order radiative corrections [8,9] or model-dependent box diagram calculation 
[10], incorporating nucleon's substructure, are not able to make a definitive conclusion on 
this discrepancy. We should note that two-photon exchange corrections are small in general 
and at a 1% level for a large class of experiments [3]. With this renewal of interest, the 
importance of theoretical descriptions of nucleon form factors in the space-like and also in 
the time-like region is emphasized. 

In principle, time-like form factors could be evaluated from the space-like equivalents by 
means of dispersion relations. The ratio between electric and magnetic proton form factors 
by using space-like and time-like data was recently analysed in the framework of dispersion 
relations [11,12]. However, all the published data in the time-like region [13-16] assumed 
Ge = Gm to hold for all and not only at threshold. As accurate data at high energy is 
lacking, and also to resolve the discrepancy between the LEAR [13] and BaBar [16] data ob- 
served close to threshold, the measurement of the proton form factors in the time-like region 
is planned at PANDA at FAIR in proton antiproton annihilation to an electron positron pair 
with unprecedented high accuracy [17]. 
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As mentioned above, the radiative corrections could be correctly evalued in the kinemat- 
ical configuration of the scattering experiment to extract the physics observables of interest. 
This paper is devoted to a theoretical investigation of the pp — > e'^e~ process including radia- 
tive corrections to lowest order of perturbation theory. Among the recent papers devoted to 
this subject, one should mention [18] where the possibility to measure the charge asymmetry 
is presented. The charge-odd part presented in the differential cross section is the origin of 
this asymmetry. It is surprising that the authors [18] do not evaluate the radiative correction 
to the cross section with an energy cut on the {e'^e~) invariant mass spectrum corresponding 
to the energy resolution of the PANDA detector. Recently, the authors of [19] reevaluate 
this correction in the laboratory frame, omitting the contribution of the hard photon as well 
as the contribution from the initial state radiation. 

The first results of full simulations with PANDA detector [17] show the precisons of the 
cross section measurement of the order of 3%-5% can be obtained. In this context, we 
need to evaluate the radiative correction due to the final state radiation, namely the photon 
emitted from the electron or the positron, as well as the radiation in the initial state at the 
proton vertex, reduced due to the large mass of the proton (antiproton) compared to the 
electron one. In this paper, the complete general radiative corrections to lowest order, to 
the pp e+e" channel are investigated , in the kinematical configuration of the planned 
experiment at FAIR with the PANDA detector. 



II. ELECTROMAGNETIC NUCLEON CURRENT OPERATOR AND BORN CROSS 
SECTION 

Let us first introduce our notations and the definition of the electromagnetic nucleon cur- 
rent operator with the magnetic and electric form factors of the nucleon. For the annihilation 

process, in the one photon exchange approximation : 

p{p~) + p{p^) e"^(g+) -t- e^{q^) 
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The corresponding Feynman diagram for this reaction is given in Fig.l. 
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FIG. 1. One photon exchange diagram for the process pp ^ e'^e 



The Born amphtude has the form 



47ra 



with 



^i(s)7. 



Fo(s) 



[7m, i] 



2) 



where the magnetic and electric form factors are related to the Dirac and Pauli form 
factors Fi and F2 : 



(1 



Ge{s)^F,{s)+tF2{s) 
Gm{s)^F^{s) + F2{s) 

-p+f , r = s/(4M2) ,131^1- AW/s 



(3) 



M and m are respectively the proton and the electron masses. 

The differential cross section in the Born approximation has the form 



da 

In 



a 



{\GM{s)\\l + cos'e) + {l-Pl)\GE{s)\'sm'e} 



(4) 



B 4:SPp 

Where 9 is the scattering angle of the positron in the center of mass system (cms) . This cross 
section was first derived by the authors of [20]. In the particular case, where the proton is 
considered as a pointlike particle, with Gm{s) = Ge{s) = 1, the formula (4) reduces to 



a 



da 



dQiB 4:s/3p 



{(2-/3|)sin2^} 



(5) 
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III. QED RADIATIVE CORRECTIONS TO FIRST ORDER 



The proton electromagnetic form factors can be extracted from the angular distribution 
of the final lepton in the Born cross section. However, this distribution is altered from 
its zeroth-order shape by radiative corrections. In practice, the distorted distribution by 
radiative effects can be written as : 



'da' 




'da' 
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idni 
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[i + S) 
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FIG. 2. Feynman diagrams for the first-order radiative correction in pp ^ e'^e . Note that the 
vacuum polarization loop (a) includes all fermions 



The set of diagrams contributing to the first order corrections are shown in Fig. 2. The 
virtual correction comes from the interference between the born diagram (Fig. 1) and the 
diagrams( (a)-(e)) of Fig. 2. The bremsstrahlung from the initial state (diagrams (f) and 
(g)) alters the effective center of mass energy and significantly changes the kinematics of 
the final lepton pair. And finally, the photon emission from the final state is represented by 
the diagrams (h) and (i). Only the bremsstrahlung corrections lead to infrared singularities. 
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These singularities are cancelled order by order by virtual corrections. We adopt the standard 
treatment of the bremsstrahlung, separating the soft photon contribution with the emitted 
photon energy up to an infrared cut-off parameter ui where the soft photon approximation 

holds, and the hard photon contribution from oo up to an experimental cut depending on the 
energy resolution of the detector. This separation is somewhat arbitrary, so we have checked 
that the total radiative correction (Virtual+soft+hard) does not depend on this infrared cut- 
off. 

The virtual and real photon corrections are achieved by the factorization of the cross section 
in Eq. (6) with 

6 = dsv + (7) 

We write down the soft and virtual correction together to remove the infrared singular- 
ities. The remaining term 5sv (soft-|-virtual) is now finite. The hard photon contribution 
5h depends, of course, on the energy resolution of the detector. The full simulations in a 
realistic configuration of the detector will allow the determination of the experimental cut 
on the maximum energy of the real photon emitted or preferably, on the invariant mass 
spectrum of the final lepton pair. 

Below we present the details of our investigation of the different contributions to the 
radiative corrections in pp — > e'^e~. 

A. Virtual and soft photon contributions 

In this section, following the previous work [18], we will suppose the proton to be the 
pointlike particle. In the range of energy and momentum transfers considered here (4M^ < 
s < 30 GeV^), the correction due to the finite size of the nucleon is found in general to be 
much smaller than the other contributions [8]. The radiative factor 5 is a ratio between the 
corrected cross section and the Born cross section where the form factors are both present. 
This could be one of the raisons, why the virtual radiative correction is pratically independant 
of the model assumption on the form factors. 
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1. Virtual correction 



The Born cross section with one-loop corrections can be written in the form 
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where s, t and u are the Mandelstam variables of the process, namely 

and is defined as the solid angle fle+ of the positron in the center of mass system. 
At lowest order in a, the vaccum correction can be written as 
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and keeping in the cross section, the correction in the same order in a, one gets 
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The real parts of different contributions through the interference with the Born amplitude 
which is real itself, are given below. 

For the vaccum polarisation contribution, the loop in diagram (a) of Fig. 2 includes all 



fermions, including the muon-loop and hadronic loop. The charged pion pair as hadron 
state is taken into account in the vacuum polarisation contribution [18]. 
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TT 
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At the lepton vertex, the expression of ^^eriex derived long time ago by the authors of 
[21]. To deal with the infrared divergent term, we consider the extra virtual photon of the 
process in Fig. 2 with a mass A. 
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Isolating explicitely the divergent term from eq. (22), one obtains : 
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At the hadron vertex, we assume the proton to be a pointlike particle and we get : 



vertex 



2Fi'^ + 
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2-/32 sin^e 2-/32 sin^^ ^ ^ 

(28) 

The expression for d^ox due to the interference between the Born amphtude and the two- 
photon exchange diagrams (d) and (e) in Fig. 2 can be found in [18] and is given in Appendix 
A. 



2. Soft photon contribution 

The soft photon contribution comes from the reaction pp — )■ 7e"'"e~ , where the energy 
of the emitted photon goes to zero. The soft photon cross section is related to the Born cross 
section using 

(io-soft = -TT^ da^ ^soft (29) 
The soft photon contribution to the cross section reads 

^R^-T^Isoit (30) 

with 

.+ ^ 2 



'soft 



k.p k.q k.q^ J 2ko 
u' = \JiJ^ — A is a virtual mass of the photon (A 0) 



4oft — 4oft + -^soft + Koi (32) 
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We adopt the 't Hooft and Veltman method [22] to evaluate these integrals, rewritten their 
final result, using our metric a la Bjorken and Drell. We give more details of this transposition 
in appendix B. 

The electron term in the contribution of the soft photon emitted at the lepton vertex can be 
written as 



-g- q- 



with 



Aw? 



2q ■ = s — 2171^ 



s - - 2 

The evaluation of the two first diagonal terms is straightforward and one gets : 
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The non-diagonal term is derived from the 't Hooft and Veltman method. We need to separate 
the finite term (finite) from the infrared singularity term (div) depending on A. 
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Now for the soft photon emitted at the hadron vertex, assuming the proton as a pointhke 
particle, /g^f^. is derived from I^^^^, replacing /3e by f3p : 



1 + ^.^,^(1 



(42) 



Here again, we separate the interference term between the soft photon contribution from 
the Icpton and hadron vertex into an infrared divergent term depending on the mass A and 
a finite term : 



^;o?(div)=4.1n^^ In— (43) 
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— 7\ In — 

M^-u M 



Sp{x) is the dilogarithm or Spence's function defined as : 



One can check that the infrared terms depending on A disappear when we sum up the 
contributions from the virtual and soft photon corrections. 

We show in table I and II, the soft and virtual corrections with final and initial state radi- 
ations. The superscript "e" or "e — p" indicate the corrections from the lepton final state 
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radiation or the total correction with final and initial state radiations, u is the infrared 
cut-off on photon energy, separating the soft and hard photon contributions. 



TABLE I. Soft and virtual corrections : s = 5.4 GeV^ , a; = 12 MeV, u/E^+ ^ 1% 



Oe^ \ S 


^vac ^vertex "vertex "_R "box W/j 


"sv "sv 


30. 

60. 

90. 
120. 
150. 


0.0103 -0.2602 -0.0147 0.0160 -0.0127 -0.0033 -0.0697 
0.0103 -0.2602 -0.0147 0.0147 -0.0127 -0.0020 -0.0365 
0.0103 -0.2602 -0.0147 0.0132 -0.0127 0.0000 0.0000 
0.0103 -0.2602 -0.0147 0.0147 -0.0127 0.0020 0.0365 
0.0103 -0.2602 -0.0147 0.0160 -0.0127 0.0033 0.0697 


-0.2646 -0.3309 
-0.2646 -0.2991 
-0.2646 -0.2641 
-0.2646 -0.2261 
-0.2646 -0.1915 


TABLE II. Soft and virtual corrections : s = 12.9 GeV^ , a; = 18 MeV, u;/E^+ « 1% 


Oe+ \ S 


<Jvac "vertex "R "vertex "R "hex 


"SV '^SV 


30. 
60. 
90. 

120. 
150. 


0.0135 -0.2991 -0.0006 0.0113 -0.0314 -0.0112 -0.1169 
0.0135 -0.2991 -0.0006 0.0095 -0.0314 -0.0059 -0.0558 
0.0135 -0.2991 -0.0006 0.0062 -0.0314 0.0000 0.0000 

0.0135 -0.2991 -0.0006 0.0095 -0.0314 0.0059 0.0558 
0.0135 -0.2991 -0.0006 0.0113 -0.0314 0.0112 0.1169 


-0.2792 -0.4162 
-0.2792 -0.3668 
-0.2792 -0.3043 

-0.2792 -0.2452 
-0.2792 -0.1824 



The radiative corrections given in the different columns are explained in the text. Let us 
make some comments on the radiative correction (soft+virtual) factor values in the two last 
columns. 

• If only the final state radiation is taken into account, it is independent of the lepton 
scattering angle in the center of mass system. This is not the case for the initial state 
radiation, due to the interference between the initial and final state radiations given by 

5r'- 
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• A remarquable feature is the asymmetry observed in the lepton angular distribution 
due to the charge-odd term. The angular distributions of the e+ and e' are different. 

• The contribution of the box diagrams (two-photon exchange) is neghgible, less than 1 
%, in agreement with the hard rescattering mechanism calculation [23]. 

In principle, the hard photon contributions in the next section do not alter these features. 
B. Hard photon emission 

We consider next the contribution from hard photon emission : 

P{p~) + p{p^)^e^{q^) + e"(g")+7(/c) 

The invariant mass W of the (e+e") system is defined as : 

W"" = (g+ + q-f = {p- + p+ - kf 

In contrast with the previous subsection, the explicit proton form factors [24] are introduced 
in our estimation of the hard photon contribution. The correction factor 5h is a ratio between 
the cross section with an extra real photon and the Born cross section and is not expected 
to be very sensitive to the explicit shape of the form factors. 

We have performed the exact calculation of 5h and the numerical results are displayed in 

Table III in terms of the total correction. 

The amplitude is written as a sum of four amplitudes : 

M^Mi + M2 + M^ + Mi (45) 

M.I, M.2, M.^, M.4 are respectively the amplitudes of the diagrams (h), (i), (f), (g) 
in Fig. 2. The expressions of different amplitudes are displayed in Appendix C. 

../^ X sr^ A^{me+,me-]mp,mp) 

M{X,me+,me-;mp,mp) = 2^ £a{kA) (46) 
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(55) 



The terms with the sum on i,j in formulas (54) and (55) are respectively the hard photon 



contributions Sjj and 5^ ^. 



TABLE III. Total radiative corrections for s = 5.4 GeV^ (two left columns) and for s = 12.9 GeV^ 
(two right columns), assuming the energy of the hard photon emission up to 100 MeV 



0e+ \ S 


6^ 




6^ 


ge-p 


30. 


-0.1154 


-0.1502 


-0.1517 


-0.2388 


60. 


-0.1154 


-0.1298 


-0.1517 


-0.1980 


90. 


-0.1154 


-0.1080 


-0.1517 


-0.1616 


120. 


-0.1154 


-0.0831 


-0.1517 


-0.1183 


150. 


-0.1154 


-0.0598 


-0.1517 


-0.0725 



The total radiative correction depends, of course, on the energy cut of the emitted photon. 
The full simulation with a reahstic configuration of the PANDA detector could determine 
the precise value of this experimental cut. 



IV. NUMERICAL RESULTS 

The numerical results strongly depend on the experimental energy cut of the emitted 
photon.This dependence is displayed in Figs. 3-4 for two values of the antiproton energy. The 
distribution in Fig. 5 allows one to construct a Monte Carlo event generator for pp — > 7e'*'e~. 
The full simulation with a realistic environnement of the PANDA detector is needed to 
determine the experimental cut on the invariant mass of the (e^e~) system. 
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FIG. 3. (Color online) Total radiative corrections factor 5 as a fonction of cos{9) in the CM frame, 
for s = 5.4 GeV^: for EJ^°-^ = 100 MeV (black line) and 50 MeV (red line). The corresponding 
dashed lines are only the final state radiation contribution 
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FIG. 4. (Color online) Total radiative corrections factor 5 as a fonction of cos{d) in the CM frame, 
for s = 12.9 GeV^: for E^™'^^' = 100 MeV (black line) and 50 MeV (red hne). The corresponding 
dashed lines are only the final state radiation contribution 
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FIG. 5. (Color online) Differential cross section (Pa/dW'^dQ,f,+ as a fonction of W'^ {W is the 
invariant mass of the (e+e") system), for s = 12.9 GeV^, E)™*^^ = 100 MeV and 6 = 30° (black 
line). The corresponding red dashed line is only the final state radiation contribution 



The corrected cross section is displayed for comparison with the Born cross section in Fig. 
6. The different curves are obtained with one model [24] for the form factors. Of course the 
normalisation of these cross sections depend on the model assumption. One can also remark 
the asymmetry of the black line due to charge-odd term when the the initial state radiation 
at the hadron vertex is included. The measurement of this asymmetry term included in the 
angular distribution seems to be a difficult task. Fig. 7 shows the ratio of the correction 
factors, reflecting the importance of the initial state radiation versus final state radiation. 
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FIG. 6. (Color online) Corrected differential cross section (i^a jd£lf,+ as a fonction of cos(0) in the 
CM frame, for s = 12.9 GeV^, E^"-^ = 100 MeV (black line). The corresponding red dashed line 
is only the final state radiation contribution. The blue dotted line is the Born cross section with a 
model assumption of the form factors [24] . 

1.2 
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cos 

FIG. 7. (Color online) Ratio between the 1 + 6^~^ and 1 + 5^ as a fonction of cos{9) in the CM 
frame, for s = 12.9 GeV^: for = 100 MeV (cyan line) and 200 MeV (red dashed line) 

The total radiative correction have an important symmetry that is worth to mention. 



18 



The values of 5 are the corrections to the angular distribution of the positron (e'*'). The 
corresponding distribution of the electron (e") is obtained by replacing (9) by (tt — 9) in 
order to respect the C-charge symmetry. We have checked numerically that the total radiative 
correction reads : 



The first one contains the charge-even terms and is the corrected Born cross section for form 
factors extraction. The second one is the charge asymmetry observable due to the odd part. 
The value of this charge symmetry term A is rather large (~ 5%) and can be measured with 
the PANDA detector. In constrast to the observable S, the charge symmetry A is more 
sensitive to the model assumption on the form factors due to the interference terms between 
the initial and final state radiations. 

V. CONCLUSIONS 

We conclude that the initial state radiation (at the hadron vertex) is not neghgible and 
can be calculated and incorporated into the total radiative correction. If the precision of the 
cross section measurement is of the order of 3-5%, the total correction obtained in this paper 
can be used to correct the Born cross section before the comparaison with the experimental 
data. We have shown that the interference between the initial and final state radiations is 
rather large. The evaluation of this correction in terms of structure functions and evolu- 
tion equations [25] must be performed with caution. A study of the angular distribution 
of the photon shows a non negligible part of the photons emitted outside the cone of angle 
9^ < m/E. 



This observation led us to define two interesting observables namely : 
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In contrast to space-like elastic electron-proton scattering, in the time-like region we can 
consider the shape of the and e" angular distributions separately to exhibit the charge 
asymmetry. The numerical result displayed in Fig. 6 shows the limit of this statement. We 

suggest the measurement of two observables namely S and A. 

We would like to emphasize that the Monte Carlo event generators have been developed 
in the course of this work. For this particular reaction, our code is more accurate than the 
PHOTOS Monte Carlo [26] where the interference between the initial and final state radiation 
is ignored. 
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Appendix A: Box diagrams contribution 

With the usual Mandelstam variables s, t and u of the process. 6 is the scattering angle 
of the positron in the center of mass system. 
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The five different terms in /(s, u) are 
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This last term contains the divergent part. Rewritting the ^\,ox in term of F^oxi one gets 
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— box 



s {u — t) 



+t + u)l,-- + -Ll- 



With 



TT 



s/3p\ ' 2 '^'^ 6 



+ 2Sp 



1 + a 



-25p 



and 



2 - box 



3 — box 



s {2t + s) 
A 

s{2u + s) 



\lI - sp 



-t 



-2Sp(- 



1-1% 



4 — boa; 



\ut- M^{s + M^)) 



-^ts -'-'us _j_ " ^ 

t U Ut ^ 



U — t 



"5 — box — 



La = In 



M2 



In 



In 



1 - 



In 



-t 



Lus = In 



M2- 



In 



M2 



(A6) 
(A7) 
(A8) 

(A9) 

(AlO) 

(All) 

(A12) 
(A13) 
(A14) 

(A15) 
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Appendix B: 't Hooft and Veltman integrals 



We need to estimate the integral of the type : 



^ f L 

J (vi ■ k) (i 



(Pk 



2 J {pi- k){pj ■ k) 2ko 
where cu' = \JliS^ — and a; = [k^max is the maximum energy of the emitted photon. 



P = VPi Q = Pj (p ~ = {r]Pi — Pj)o and pj^ with the same sign 



v'^Pi - '^VPi-Pj +P'j = 



p2 _ q2 

2£ 



Cij {div) = 2tt-^ In ^ In ^ 



£jj (finite) = 27r 



ve 



l|^2 P0-|Pl 1^^2g0-|g| 



4 po+ \p\ 4 qo + \q\ 



+27r 



Sp(l-P^)-Sp{l-'-^) 

sp(i-p^ysp{i-'-^) 



Appendix C: Hard photon contribution 



Defining 



iH-l^^) F^o{e)^l F2o{e)^K, 



q = p +p^ — k = q — k 
22 



Mi{X,me+,me-;mp,mp) = — — — ^ — ^ ^aV^^^) 



(C3) 



(me+ , me- ; rup, rup) = ^ el- Vp (q) Up 



4+ + m)Yve+ (C4) 



Di = {k + q-^f -m^ ^2k-q 



-,+ 



(C5) 



M2{X, me+ , me- ; mp, nip) 



A^{me+,me^;mp,mp) 



(C6) 



A'^{me+,me-;mp,mp) = —Cpcl- VpT%j^Jq)up Ue- + i~ + m) j^Ve+ (C7) 



i:)2 = {k + q~y -m^ ^2k-q 



(C8) 



3 ( A, me+ , TOe- ]mp,mp) 



A^{me+,me-;mp,mp) 



e:{k,X) (C9) 



^3(me+,me-;mp,mp) 



(CIO) 



^3 = {p- -ky-M^ = -2k-p 



(Cll) 



. , /, s Al[mg+,me-;mp,mp) 
M4{\me+,me-;mp,mp) = — — e^{k,X) 

JJa 



(C12) 



Al (me+ , rUe- ; rup, mp) = 



V2 S^e- 



Ue-lfjLVe+ 



(C13) 
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^4 = (p+ -kf-M^ = -2k- 



(C14) 
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